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一 种 尺度 空间 特征 区 域 的 强 鲁 棒 性 水 印 算法 
齐 向 明 ， 徐 ” 召 ， 李 ” 表 ， 侯 明 帮 


(辽宁 工程 技术 大 学 软件 学 院 , 辽宁 戎 芦 岛 125100) 


摘 要 : 针对 数字 水 印 算法 中 特征 区 域 选取 不 足以 反映 图 像 重要 信息 ， 导 致 便 棒 性 减弱 的 问题 ， 提 出 一 种 尺度 空间 特 
征 区 域 的 强 鲁 棒 性 水 印 算法 。 通 过 尺度 空间 特征 点 检测 ， 提 取 靠 近 载 体 图 像 重 心 且 互 不 重 登 的 特征 区 域 ， 合 成 特征 区 
域 和 矩阵 ， 用 变换 域 水 印 算 法 嵌入 水 印 ， 对 可 能 受到 攻击 的 含水 印 图 像 提 取 特 征 区 域 ， 合 成 其 特征 区 域 矩 阵 ， 用 谋 入 水 
印 算 法 逆 过 程 提取 水 印 。 实 验 表明 ， 该 算法 不 仅 对 一 系列 攻击 都 具有 强 鲁 棒 性 ， 且 不 可 见 性 良好 。 
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Strong robust watermark algorithm in scale space signature area 


Qi Xiangming, Xu Man, Li Yue, Hou Mingjun 
(College of Software Liaoning Technical University, Huludao Liaoning 125105, China) 


Abstract: Area selection based on the digital watermarking algorithm is not enough to reflect the image ofimportant information, 

lead to decreased robustness problems, this paper put forward a measure space characteristics of the area the strong robustness 

of watermarking algorithm. Through scale space feature point detection and extraction characteristics near the image center of 

gravity and non-overlapping regions, regional matrix synthesis characteristics, with the transform domain watermarking 

algorithm embedding watermarking, the area may contain the watermark image feature extraction under attack, the 

characteristics of regional matrix synthesis, it used inverse process to extract the watermark embedding watermarking algorithm. 

Experiments show that this algorithm is not only robust to a series of attacks, but also has good invisibility. 

Key words: digital watermarking; scale space; characteristic region; transform domain; strong robustness 

于 和 提取 ， 特 别 是 辅助 确定 水 印 衣 入 的 准确 位 置 ， 有 助 于 增强 水 

0 印 算法 的 鲁 棒 性 。 文 献 [8,9] 对 采用 Harris 算 子 检测 特征 点 ， 实 

网 络 的 开放 性 与 共享 性 给 人 们 生产 、 工作 和 生活 带 来 了 极 ” 现 对 图 像 重 要 特征 的 提取 ; 文献 [10] 在 小 波 变 换 域 提取 的 Harris 
大 地 便利 ， 随 之 也 产生 了 媒体 信息 的 安全 问题 。 数 字 水 印 技术 ”特征 点 嵌入 水 印信 息 ， 在 抗 几何 攻击 上 体现 了 较 好 的 鲁 棒 性 ; 
作为 知识 产权 保护 和 认证 的 一 种 有 效 方式 ， 得 到 了 广泛 应 用 。 文献 [11] 提 出 了 基于 直方 图 均衡 化 和 Harris 角 点 检测 的 数字 水 
水 印 算法 的 关键 是 在 良好 不 可 见 性 的 前 提 下 ,增强 水 印 鲁 棒 性 。 ” 印 算 法 ,可 提取 数量 较为 可 观 的 特征 点 ;文献 [12,13] 利 用 Harris- 
强 鲁 棒 性 水 印 算法 成 为 国内 外 专家 和 学 者 研究 的 热点 对 象 。 Affine 提取 出 图 像 的 仿 射 不 变 特征 点 ， 特 征 点 稳定 性 较 好 ; 

变换 域 较 空域 更 容易 设计 出 强 鲁 棒 性 的 水 印 算法 中。 基于 ”Mikolajczky 等 人 0 将 Harris 特征 点 检测 算 子 与 尺度 选择 相 结 
变换 域 的 水 印 算法 主要 通过 修改 变换 的 系数 值 存储 水 印 外 。 文 ” 合 , 提出 了 Harris-Laplace 特征 点 检测 算 子 , 既 能 检测 到 较 大 数 
献 [3~5] 将 载体 图 像 进行 单 一 变换 ,算法 只 能 在 对 抗 某 一 种 攻击 。 量 的 特征 点 ， 且 包含 尺度 特征 。 综 上 所 述 的 各 种 基于 图 像 特 征 
上 表现 较 好 的 鲁 棒 性 ， 且 其 不 可 见 性 较 差 ， 文献 四 将 载体 图 像 往 域 水 印 算法 ， 虽 在 某 方面 得 到 改进 ， 但 没有 兼顾 解决 以 下 几 
进行 DCT 变换 和 SVD， 一 定 程 度 均 衡 了 水 印 算法 的 不 可 见 性 ”点 问题 ，a) 所 提取 图 像 特征 点 数量 较 少 ， 不 能 较 好 地 反映 图 像 
和 和 鲁 棒 性 ， 但 鲁 棒 性 仍然 较 差 ， 文献 [7] 将 载体 图 像 先 作 DWT ”的 某 些 特征 : b) 特 征 点 检测 方法 存在 很 多 元 余 点 ， 增 加 嵌入 水 
变换 再 作 SVD， 得 到 载体 图 像 奇异 值 矩阵 ， 将 置 乱 后 的 水 印 作 。” 印 时 间 且 错误 率 较 高 ，c) 特 征 点 稳定 性 差 且 分 布 不 均匀 ， 影 响 
相同 处 理 得 到 水 印 奇异 值 矩阵 ， 把 后 者 加 权 嵌 入 到 前 者 较 大 系 ”数字 水 印 抗 攻 击 能 力 ; dd) 不 能 结合 图 像 内 容 确 定局 部 特征 区 域 
数 中 ， 载 体 图 像 奇 异 值 不 发 生 剧 烈 改变 ， 但 水 印 嵌 入 不 完全 ， 尺寸 和 位 置 ， 减 弱 了 水 印 算法 的 鲁 棒 性 。 
鲁 棒 性 不 理想 。 通 过 对 图 像 本 身 的 某 些 特征 来 辅助 水 印 的 嵌入 基于 以 上 分 析 , 本 文采 用 改变 Harris-Laplace 算法 中 Harris 
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角 点 的 尺度 空间 来 获得 特征 点 ， 根 据 水 印 图 像 大 小 确定 特征 区 
域 大 小 ,选取 载体 图 像 中 靠近 图 像 重心 且 互 不 重 登 的 特征 


x 域 ， 


着 图 像 的 能 量 特性 ， 当 图 像 受到 攻击 时 , 变换 后 DCT 的 低频 系 
数 和 DWT 后 的 低频 子 带 系数 具有 相同 的 变化 趋势 ， 即 能 同时 
变 大 或 变 小 19。 因此, 利用 关系 双方 DCT 和 DWT 系数 上 述 特 


合成 特征 区 域 矩阵 ， 对 其 依次 进行 DWT、DCT 和 SVD 变换 ， 
完成 水 印信 息 甬 入 。 对 可 能 受到 攻击 的 含水 印 图 像 进行 特征 区 
或 的 合成 ， 再 作 水 印信 息 舱 入 的 逆 过 程 ， 提 取水 印 。 该 算法 的 
标 是 增强 各 类 攻击 鲁 棒 性 ， 同 时 不 影响 水 印 的 不 可 见 性 。 


1 ”基础 理论 


1.1 尺度 空间 特征 点 检测 
尺度 空间 理论 的 基本 思想 是 : 在 图 像 信 息 处 理 过 程 中 引入 
一 个 被 视 为 尺度 的 参量 ， 通 过 不 断 变化 参量 获得 多 尺度 下 的 空 
间 序 列 ， 对 这 些 序列 进行 尺度 空间 主 轮 廓 的 提取 ， 并 以 该 主 轮 
廓 作为 特征 向 量 ， 实 现 对 边缘 、 角 点 检测 和 不 同 分 辨 率 上 的 特 

。Harris-Laplace 特征 点 检测 算 子 计算 特征 点 的 尺度 信息 ， 
尺度 信息 通常 是 未 知 的 ， 一 般 而 言 ， 较 大 尺度 能 较 可 靠 地 消除 
误 检 ， 但 定位 不 够 准确 ; 反之 ， 较 小 尺度 能 准确 定位 真正 的 特 
征 点 ， 而 误 检 率 较 高 。 因 此 ， 通 过 引入 尺度 空间 理论 ， 在 多 尺 
度 空 间 上 寻求 最 大 响应 点 ， 即 该 点 在 正 交 方向 上 的 梯度 变化 显 


著 ， 将 其 作为 候选 特征 点 ， 求 得 候选 特征 点 集 { 严 }。 


确定 候选 特征 点 集 基 础 上 ， 计 算 每 一 特征 点 的 特征 尺 
(characteristic scale), 它 反映 局 部 图 像 特 征 与 操作 算 子 间 的 最 大 
相似 度 ， 选 用 LOG 函数 计算 特征 尺度 。LOG 算 子 如 下 : 
OG(%,y,60) ,0 GH y,60),) 

Ox” Oy” 


洱 


LOG(x, y,6,) = 62 


其 中 : 数字 图 像 记 为 1 对 于 给 定 的 尺度 空间 6 =1.4"6，， 


59 =0.76 (n=1,2,.…….,15) 及 闵 值 1 =1000， 对 于 


{ 及 } 中 的 每 个 候选 特征 点 ， 和 迭代 确定 最 终 的 图 像 特 征 点 与 特 


征 尺度 。 实 现 方法 如 下 : 
a) 检验 LOG 算 子 在 预选 的 图 像 特 征 点 处 ， 能 否 在 尺度 搜 
索 范 围 内 获得 局 部 极 值 ， 如 果 得 不 到 极 值 ， 该 点 舍弃 。 尺 度 搜 


其 中 f= 二 0.7,...,1.4. 


索 范围 限定 为 :6% =1504) ， 


b) 将 LOG 算 子 取 极 值 的 特征 点 及 加 入 候选 特征 点 集合 


{P.}. 


复 以 上 两 步 , 至 所 有 特征 点 处 理 完毕 为 止 .集合 { 己 } 


le) 

b= 
Im 
Mm 


就 是 所 求 Harris-Laplace 特征 点 。 
1.2 ， 强 鲁 棒 性 水 印 嵌 入 和 提取 
小 波 变换 后 的 低频 子 带 和 余弦 变换 后 的 低频 系数 都 表征 


征 进行 的 水 印 算法 具有 一 定 的 鲁 棒 性 。 
从 图 像 处 理 的 角度 来 看 ,奇异 值 分 解 特性 是 ,一 幅 图 像 的 
奇异 值 具 有 相当 好 的 稳定 性 ， 即 当 图 像 受到 轻微 扰动 时 ， 其 奇 
异 值 不 会 发 生 剧烈 的 改变 ， 奇 异 值 能 够 表现 出 图 像 内 在 的 代数 


特性 09。 大 小 为 HXn 的 年 阵 A E Ce 的 奇异 值 分 解 为 


A=UDV” = > A4uv” 
El 


其 中 : Ue C”™ ,Ve C”™Y 是 两 个 西 矩 阵 , 即 UT?U = 瓦 ， 


VV = 万 ;DD 表示 元 素 为 非 负 数 的 对 角 阵 , 7 < min{7m,n} 


阵 A 的 秩 ， 岂 = VU 称 为 矩阵 A 的 奇异 值 ，J4; 


Pa 
> 
[OIL 


(1 = 1,2,...,7) 称 为 矩阵 A 的 左 奇异 向 量 ， 对 应 于 矩阵 AA7 


的 特征 向 量 ; 


V, (1 = ],2,.….,7) 称 为 矩阵 A 的 右 奇异 向 量 ， 


对 应 于 矩阵 A7 4 的 特征 向 量 。 
像 矩 阵 的 奇异 值 与 图 像 内 在 的 特征 存在 相关 性 5 , 具 
较 强 的 稳定 性 。 当 图 像 遭 受到 较 小 的 攻击 时 ， 其 奇异 值 不 会 发 
生 较 大 的 改变 , 且 图 像 的 奇异 值 对 图 像 失真 具有 不 变性 。 因 此， 
将 SVD 与 DWT 和 DCT 相 结合 ,算法 的 时 间 复 杂 度 大 大 降低 ， 
水 印 具 有 更 好 的 鲁 棒 性 和 不 可 见 性 。 


pa 


2 ”算法 实现 


2.1 图 像 预 处 理 

2.1.1 水 印 图 像 加 密 
为 消除 像素 空间 的 相关 性 , 需要 对 水 印 图 像 进 行 加 密 处 理 。 

采用 Logistic 映射 混沌 加 密 作为 处 理 方法 。Logistic 映射 定义 如 


Kin = Xe XAX(—X) 


其 中 ， 4 E [0,4] 被 称 为 Logistic 参数 ; NX E [0,1]， X0 即 


为 密 钥 key。 在 当 3.5699456< 4 < 4 时 ，Logistic 映射 处 
于 混沌 状态 。 处 于 混沌 状态 的 水 印 图 像 被 随机 而 均匀 地 置 乱 ， 
产生 了 一 个 非 周 期 、 不 收敛 的 ， 并 对 初始 条 件 敏 感 的 序列 ， 像 
素 间 相关 性 被 消除 。 
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2.1.2 载体 图 像 特征 区 域 选 取 

通过 Harris-Laplace 算 子 提取 图 像 的 特征 点 。 载 体 图 像 
Lena、Barbara、Baboon、door 和 flower 进行 Harris-Laplace 特 
征 点 检测 的 结果 见 图 1。 可 以 看 出 ， 特 征 点 的 分 布 区 域 主要 在 
复杂 纹理 和 灰 度 变化 明显 的 位 置 ， 不 是 所 有 的 特征 点 都 适用 于 
谍 入 水 印 。 本 算法 以 特征 点 为 中 心 ， 选 取 靠 近 载 体重 心 的 特征 
区 域 ， 特 征 区 域 可 能 存在 重 炙 。 为 避免 水 印信 息 互相 影响 ， 需 


要 筷 选 出 部 分 不 相 重 准 的 特征 区 域 作为 最 佳 的 肉 入 水 印 的 位 置 。 


筛选 的 步骤 如 下 所 述 


图 1 Harris-Laplace 特征 点 检测 结果 


a) 载体 图 像 1I 通 过 尺度 空间 的 Harris-Laplace 算 子 得 到 特 


征 点 集合 { 已 }; 


b) 根据 水 印 图 像 的 大 小 ， 选 择 特征 


区 域 的 大 小 。 以 { P, } 


中 特征 点 为 中 心 取 抢 形 特征 区 域 ; 


内 
IN 
ua 


c) 从 { 己 } 中 删 掉 特征 点 对 应 的 矩形 特征 区 域 存在 有 


部 分 的 特征 点 ， 并 保存 为 { 已 } 


d) 标记 载体 图 像 的 重心 位 置 { 已 } 中 每 


[1(X11, yi)， 集 合 


个 特征 点 1 的 距离 记 做 集 { Dp }， 对 { Dp } 从 小 到 大 排序 。 


e) 根据 水 印 图 像 的 大 小 , 选择 特征 区 域 的 个 数 , 依次 按照 


{Dp } 中 的 顺序 取 对 应 的 特征 点 及 特征 区 域 作为 水 印 嵌入 子 


块 ， 保 证 了 特征 区 
2.2 水 印 嵌入 
设 原 始 载 体 图 像 I 的 大 小 为 M x M ， 置 乱 后 的 水 印 
像 为 WW， 其 大 小 为 Nx N 。 水 印 杠 入 的 具体 流程 如 下 : 
a) 将 特征 区 域 合成 特征 区 域 矩 阵 AA , 作 离 散 小 波 变换 ， 


域 靠 近 载 体 图像 的 重心 


且 互 不 重奏 。 


| 


对 低频 子 带 LL 离散 余弦 变换 得 到 DCT 系数 矩阵 上 Ly， 奇异 


值 分 解 后 ( 即 LL; = UDV”) 形 成 正 交 


矩阵 U 、W 和 对 角钢 


[GT 


Ri 全 优 Fi, 


阵 8S=diag(1 ,1 )。 


b) 通 过 式 Temp = 9 十 Qf *W 修改 载体 图 像 对 角 阵 S ， 


其 中 ，af 为 嵌入 强度 , 本 算法 中 取 Qf = 47 . 对 矩阵 Temp 


作 奇 异 值 分 解 ， 得 到 正 交 和 矩阵 Ul、V1l 和 对 角 和 矩阵 


S1= diag(1, 力 ,.….,Av)， 用 Sl 替换 S 进行 逆 SVD 变 


换 , 得 到 杠 入 水 印 后 的 DCT 系数 矩阵 , 对 其 进行 反 离散 余弦 变 
换 ， 将 其 低频 子 带 和 其 他 中 高 频 方 向 子 带 作 逆 小 波 变换 ， 得 到 
含有 水 印 的 特征 区 域 矩 阵 AA 。 
c) 对 含有 水 印 的 特征 区 域 矩阵 AA 进行 分 解 ， 依 次 蔡 换 
原 图 像 中 的 特征 区 域 ， 得 到 含有 水 印 的 图 像 了 。 
2.3 水 印 提取 
水 印 提取 前 , 对 可 能 受到 攻击 的 含 
成 其 特征 区 域 矩 阵 ， 对 该 矩阵 使 ) 
息 。 水 印 提取 的 县 体 流程 如 下 : 
a) 提 取 可 能 受到 攻击 的 含水 印 图 像 了 的 特征 区 域 , 对 特征 
区 域 进行 合成 运算 ， 构 成 可 能 受到 攻击 嵌入 水 印 特征 区 域 矩阵 
BB. 
b) 对 BB 进行 离散 小 波 变 换 ， 对 低频 子 带 LL1 进行 离散 
余弦 变换 ， 将 得 到 的 DCT 系数 矩阵 记 为 LLL ， 对 LLL 作 奇 
异 值 分 解 ， 得 到 正 交 和 矩阵 U2 或 V2 以 及 对 角 和 矩阵 


水 印 图 像 提取 特征 区 域 ， 
j 拒 入 水 印 逆 过 程 提 取水 印 


zm > 


S2= diag(4,h,..., Ny): 


co) 用 5S2 替换 Sl 进行 逆 SVD 变换 ， 得 到 中 间 和 矩阵 SN ， 


SN—S 


通过 式 WN = 提取 出 水 印信 息 WN 。 


dd) 根据 多 = 二 4、Xo = 0.2345 通 过 Logistic 映射 混沌 解 


密 将 提取 出 的 水 印信 
图 像 WNN 。 


息 WAN 进行 解密 运算 , 得 到 解密 后 的 水 印 


3 ”仿真 实验 及 分 析 


3.1 实验 环境 及 参数 说 明 

仿真 实验 采用 MATLAB R2015b 作为 实验 平台 ， 选 取 
Lena、Barbara 和 Baboon 三 幅 标准 灰 度 图 像 作 为 载体 图 像 ， 其 
大 小 为 384x384。 选 取 大 小 为 32x32 的 二 值 图 像 “ 辽 宁 工 大 ” 


为 水 印 图像 。 当 XeE[0,1]， 多 越 接 近 于 4 时 ，Logistic 映射 


工作 越 处 于 混沌 状态 , 因此 取 Logistic 映射 混沌 加 密 中 Logistic 


参数 4 = 4 , 密 钥 随机 取 值 为 Y0 = 0.2345 . 图 2(a) 是 32x32 
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的 二 值 水 印 图 象 ， 对 该 水 印 进 行 Logistic 映射 混沌 加 密 ; 


水 印 图 像 经 过 混沌 加 密 的 加 密 图 像 。 


(a) 水 印 图 像 (b) 加 密 图 像 

图 2 水 印 图像 经 过 Logistic 映射 混沌 加 
采用 峰值 信 品 比 (PNSR) 来 评价 原始 图 

间 的 质量 差别 ， 同 时 采用 归 一 化 互相 关系 数 (NC) 来 订 


密 后 的 加 密 图 
像 和 含 


像 


(b) 是 


水 印 图 像 之 
F 价 原始 水 


印 与 提取 解密 的 水 印 的 相似 程度 , NC 值 越 大 , 表示 两 者 相似 程 


度 越 大 ， 则 水 印 的 鲁 棱 性 越 强 。 
3.2 ”和 鲁 棒 性 测试 


为 了 检测 本 算法 具有 较 好 的 鲁 棒 性 , 对 含 7 
进行 噪声 攻击 、 旋 转 攻击 、 缩 放 攻击 等 多 种 类 型 


水 印 的 


且 在 水 印 提取 后 对 水 印 进行 Logistic 映射 混沌 解密 。 


3.2.1 嗓 声 攻击 


向 五 幅 含 水 印 的 图 像 分 别 添加 乘积 


盐 噪声 和 高 斯 白 噪 声 。 表 1 列 出 了 五 


必 


生 噪 声 、 高 


幅 图 像 添加 噪声 攻击 后 提 


图 像 分 别 


4 攻击 实验 ， 


[噪声 、 椒 


取出 的 水 印 NC 值 。 尺 度 空间 的 特征 区 域 均 靠 近 载 体 图 像 的 重 


心 位 置 ， 反 映 载 体 图 像 的 较 多 信息 ， 
小 波 变换 ， 并 对 其 低频 子 带 进行 离散 
数 和 矩阵 进行 奇异 值 分 解 ， 而 加 噪 是 对 
乎 不 会 影响 到 在 特征 区 域 低频 中 肉 入 
中 内 入 的 水 印信 息 产生 一 定 的 干扰 。 


提高 水 印 抗 乘积 性 噪声 、 高 斯 噪声 、 椒 盐 噪 声 和 离 


通过 对 特征 
余弦 变换 ， 得 至 
图 像 高 频 信 息 的 


区 域 进行 离散 
| 的 DCT 系 
操作 ， 几 


的 水 印信 息 ， 只 会 对 高 频 
忆 此 ， 本 文 算 法 可 


以 有 效 


斯 


噪声 攻 


总 


击 的 能 力 。 在 嵌入 水 印 前 对 图 像 进行 Logistic 映射 混沌 加 密 ， 


进一步 分 散 水 印信 息 ， 使 其 能 提取 


有 较 大 NC 


值 的 水 印 。 图 
3 展示 了 各 种 噪声 攻击 下 的 含水 印 图 像 及 其 所 提取 出 的 水 印 图 


像 。 图 中 的 水 印 图 像 清晰 可 辨 ， 说 明 本 算法 可 以 有 效 抵抗 噪声 


攻击 ， 具 有 强 鲁 棒 性 。 
表 1 含水 印 图 像 添加 噪声 攻击 后 


提取 出 的 水 印 NC 值 


齐 向 明 ， 等 


(c2) 


(c3) 


(c4) 


(c5) 
(c) 乘 积 性 噪 


一 种 尺度 空 hn 全 优 


辽宁 
工大 


(d1) 


辽宁 
二 六 


(d2) 


辽宁 
工大 


(d3) 


辽宁 
二 大 


(d4) 


辽宁 
二 六 


(d5) 


声 (0.02) 


3.2.2 旋转 攻击 


对 五 幅 含 水 印 的 图 像 分 别 进行 向 左 或 向 
。 表 2 列 出 了 旋转 攻击 后 提取 出 的 水 印 NC 值 。 对 于 


图 3 噪声 攻击 实验 图 像 


Fen 


(qd) 提取 解密 的 水 印 


右 多 种 幅度 的 旋 


F 不 同 角 


攻击 方式 参数 Lena Barbara Baboon door flower 
乘积 性 噪声 。 0.01 .0000 0.9972 1.0000 1.0000 1.0000 
0.015 1.0000 1.0000 0.9972 0.9944 1.0000 
0.02 0.9944 0.9972 1.0000 0.9916 0.9911 
高 斯 噪声 0.01 0.9944 0.9944 0.9917 0.9916 0.9930 
0.005 1.0000 0.9972 0.9972 1.0000 0.9986 
椒盐 噪声 0.01 1.0000 “0.9972 1.0000 1.0000 0.9993 
0.015 0.9972 0.9972 0.9972 1.0000 0.9979 
高 斯 白 噪 10 1.0000 1.0000 1.0000 1.0000 1.0000 
15 1.0000 1.0000 1.0000 1.0000 1.0000 


0.9944 


0.9945 0.9916 0.9944 0.9937 


党 炮 


、 不 同方 向 的 旋转 攻击 ， 提 取出 的 水 印 NC 值 均 在 0.9700 以 
。 图 4 是 旋转 角度 为 30"、 方 向 为 向 左 时 旋转 攻击 攻击 下 含 


水 印 图 像 及 其 所 提取 出 的 水 印 图 像 。 由 于 图 中 水 印 图 像 足够 清 
晰 ， 表 明 本 算法 对 旋转 攻击 可 以 有 效 抵抗 ， 鲁 棒 性 强 。 


表 2 含水 印 图 像 旋转 攻击 后 提取 出 的 水 印 NC 值 


尺 y 合 作 期 和 


hi 
齐 向 明 ， 等 : 一 种 尺度 空间 特征 区 域 的 强 前 棒 性 
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法 
参数 Baboo 
问 左 0.9945 0.9861 0.9944 0.9888 0.991 条 宁 
。 下 志 
问 右 。 0.9974 0.989 1.0000 0.9972 0.9954 
10 
问 左 0.9873 0.9725 0.9945 0.9862 0.9844 (el) GD) 
20 (e) 向 左旋 转 30 (提取 解 密 的 水 印 
句 右 。 0.9903 0.9863 0.9945 0.9861 0.9893 图 4 ”旋转 攻击 实验 图 像 
20 3.2.3 压缩 攻击 
句 左 0.9861 0.989 0.9945 0.9751 0.9862 为 检测 算法 对 不 同 强 度 JPEG 压缩 攻击 抵抗 能 力 ， 对 五 幅 
30 含水 印 的 图 像 分 别 进 行 一 系列 压缩 攻击 。 表 3 列 出 的 是 图 像 经 
句 右 。 0.9889 0.9945 0.9917 0.9806 压缩 攻击 后 提取 出 的 水 印 NC 值 。 对 于 不 同 压缩 强度 的 攻击 ， 
30 WW 提取 出 的 水 印 NC 值 均 在 0.980 0 以 上 。 图 5 是 压缩 质量 因子 为 
30% 时 的 压缩 攻击 下 含水 印 图 像 及 其 所 提取 出 的 水 印 图 像 。 图 
中 的 水 印 图 像 清晰 可 辨 ， 表 明 本 算法 可 以 有 效 抵 抗 噪声 攻击 ， 
具有 强 鲁 棒 性 。 
订 宁 表 3 含水 印 图 像 JEPG 压缩 攻击 后 提取 出 的 水 印 NC 值 
工大 攻击 方式 参数 Lena Barara 0 door flower 
n 
0.991 
10 0.9970 1.0000 0.9872 0.9809 
(el) (1) 
0.994 
20 0.9972 1.0000 0.9972 0.9836 有， 
pb a 0.996 
TT 压缩 攻击 30 0.9972 1.0000 1.0000 0.9889 ， 
二 六 0.997 
50 1.0000 1.0000 1.0000 “0.9889 
(el) (GD 70 1.0000 1.0000 10000 1.0000 ee 
1.000 


100 1.0000 1.0000 1.0000 1.0000 


(el) (f1) 
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(g1) 


(gl) 


(BJEPG 压缩 (30) 


图 5 


3.2.4 其 他 类 型 攻击 


旋转 攻击 实 


(h) 


(h) 


J 本 
工大 


二 | 
1 


(h) 提 取 解 密 的 水 印 
验 图 像 


齐 向 明 ， 等 : 一 种 尺度 空间 


品 比 PSNR 值 ，(k) 为 提取 的 水 印信 息 ; 
息 ， 分 别 是 两 组 不 同 尺寸 大 小 的 水 印 ， 并 标 


值 。 


仿真 实验 表明 ， 本 文 提出 的 水 印 算法 
匀 衡 不 可 见 性 。 从 视 


较 好 地 ] 
图 像 没有 明显 变化 ， 


Na 全 但 性 期 刊 。 


(1) 为 解密 后 的 水 印信 
Fk 注 计算 得 到 的 NC 


具有 强 鲁 棒 性 


ba 同时 


觉 效 果 来 看 ， 


辨 还原 度 很 高 。 


上 ， 具 有 较 好 的 感官 效果 ， 


水 印 的 不 可 见 性 


含水 印 的 图 像 与 原始 
图 像 质量 良好 ， 并 且 提取 出 的 水 印 清晰 可 
算法 得 到 的 峰值 信 噪 比 PSNR 都 达到 45 以 
E 较 好 ， 提 取出 来 的 


水 印 图 像 NC 值 都 达到 1.000 0 以 上 ， 水 印 具 有 强 鲁 棒 性 。 实 验 


结果 验证 了 本 文 算法 在 水 印 良好 不 可 见 性 


鲁 棒 性 ， 具 有 较 好 的 实用 性 和 合理 性 。 


的 前 提 下 ， 实 现 了 强 


对 含水 印 的 图 像 分 别 进 行 图 像 变 上 暗 攻击 、 增 加 对 比 度 、 高 
斯 低 通 滤波 攻击 和 缩放 攻击 。 表 4 是 图 像 在 受到 不 同 攻击 后 所 
提取 水 印 的 NC 值 。 随 着 攻击 强度 的 增 大 , 对 提取 出 的 水 印 NC 
值 有 一 定 的 影响 ， 但 仍 能 提取 出 具有 较 大 NC 值 的 水 印 。 本 算 

法 可 以 有 效 抵抗 各 种 类 型 的 攻击 ， 具 有 较 强 的 鲁 棒 性 。 

表 4 含水 印 图 像 受 攻击 后 提取 出 水 印 NC 值 
参数 Lena Barara Baboo door flower 
攻击 方式 
n 
(0,0.98) 0.9972 1 1 0.9991 ”0.9990 
图 像 变 瞳 (0,0.95) 0.9916 0.9916 0.9944 0.9925 0.9925 
(0,0.92) 0.9861 0.9888 0.9888 0.9879 0.9879 
高 斯 低 通 滤波 3x3 0.9945 1 0.9972 0.9972 0.9972 
增加 对 比 度 (0.1,0.9 0.9833 0.9917 0.9833 0.9861 0.9861 
) 

缩放 攻击 1/2 0.9848 0.9770 0.9931 0.9890 0.9889 
1.2 0.9975 0.9832 0.9956 0.9966 0.9965 


3.3 不 可 见 性 性 测试 


在 载体 图 像 上 分 别 


像 。 图 6G) 为 标准 灰色 
区 入 水 印 后 的 载体 


flower; 0) 为 


嵌入 不 同 大 小 水 印信 息 得 到 的 含水 印 图 
图 像 Lena、Barbara、Baboon 、 


图 像 ， 并 标 有 该 


door 和 
图 像 的 峰值 信 


PNSR=45.0579dB 
i 
二 宁 
工大 


NC=1.0000 


PNSR=45.2302dB 
1 一 一 
1 了 本 
工本 


NC=1.0000 
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录用 稿 齐 向 明 ， 等 : 一 种 尺度 空间 hoa 人 但 
高 斯 低 通 滤波 3X3 0.82 0.9734 ”0.9420 0.9945 
增加 对 比 度 (0.1,0.9) 0.9634 0.9521 0.9872 0.9833 
door PNSR=46.5898dB 和 
工大 
NC=1.0000 
杂 委 性 噪声 高 亲信 产 向 抹 ET ee a 图 像 变 暗 缩放 攻击 ”高 斯 低 遂 增强 对 比 度 
图 7 Lena 图 像 不 同 算法 性 能 比较 折线 图 
为 验证 本 文 水 印 算法 同时 良好 的 不 可 见 性 ,选择 含水 
印 的 Lena 图 像 的 PSNR 值 与 文献 [18~20] 进 行 对 比 , 本 文 算法 、 
a 文献 [18~20] 得 到 的 PSNR 值 分 别 为 47.302 6、40.290、45.700、 
订 宁 42.47， 且 本 文 算法 得 到 的 含水 印 图 像 的 PSNR 值 都 能 达到 45 
工 未 以 上 ， 可 见 本 文 算法 得 到 的 含水 印 图 像 具 有 更 好 的 不 可 见 性 。 
与 文献 [18~20] 对 比 ， 本 文 算法 实现 更 精准 的 水 印 提取 解密 ， 对 
Ee 于 各 种 攻击 的 抵抗 能 力 都 保持 在 一 个 较 高 的 水 平 ， 可 见 本 文 算 
法 能 够 较 好 地 增强 水 印 算法 的 鲁 棒 性 ， 提 高 水 印 算法 抵抗 多 种 
人 攻击 的 能 力 ， 同 时 均衡 水 印 的 不 可 见 性 。 
(9 提取 出 水 印 。 0) 解密 后 水 印 
图 6 载体 图 像 、 嵌 入 水 印 图 象 、 提 取出 水 印 以 及 解密 后 水 印 图 像 4 ”结束 语 
3.4 对 比 验 证 
为 验证 本 算法 具有 较 强 的 鲁 棒 性 ， 选 择 Lena 图 像 为 载体 rr Laplaee FeAl A nT 
图 像 的 实验 结果 与 文献 [18~20] 算 法 进行 对 比 ， 其 NC 值 对 比 实 征 点 ， 根 据 水 印 图 像 的 大 小 确定 特征 区 域 大 小 ， 选 择 靠近 载体 
验 结果 如 表 5 所 示 。Lena 图 像 不 同 算法 性 能 比较 折线 图 如 图 7 重心 且 互 不 重 登 的 特征 区 域 作 为 待 嵌入 水 印 区 域 ， 实 现 水 印 的 
所 示 。 从 表 5 及 图 7 可 以 看 出 ， 与 文献 [18] 相 比 ， 本 文 在 受到 重要 特征 区 域 贡 入 。 在 水 印 能 入 过 程 中 ， 对 特征 区 域 矩 阵 作 离 
噪声 攻击 、 旋 转 攻 击 、 压 缩 攻击 以 及 其 他 多 种 类 型 攻击 后 所 提 ”” 散 小 波 变换 和 离散 余弦 变换 ， 使 水 印 可 以 嵌入 到 特征 区 域 的 低 
取出 的 水 印 NC 值得 到 了 显著 提升 ， 说 明 其 在 对 抗 攻击 方面 的 。” 诺 中 ， 水 印信 息 的 嵌入 和 提取 同步 ， 结 合 奇异 值 分解， 提高 水 
优势 突出 ， 鲁 棒 性 能 明显 优 于 文献 [19,20] 的 算法 。 印 算法 的 稳定 性 。 在 水 印 提取 前 ， 对 相应 特征 区 域 进行 水 印信 
息 提取 。 仿真 实验 表明 ，Harris-Laplace 特征 点 检测 算 子 与 变换 
de 域 水 印 算法 结合 ， 水 印信 息 被 嵌入 到 载体 图 像 的 重要 特征 区 域 
a 中 ， 兼 顾 解决 了 原 有 水 印 算法 在 特征 点 检测 和 特征 区 域 选 取 过 
U9 09 P90 和 1 流程 中 存在 的 问题 。 基 于 尺度 空间 特征 区 域 的 强 鲁 棒 性 水 印 算法 
WA 能 够 有 效 抵抗 多 种 类 型 的 攻击 ， 实 现 了 水 印 的 强 鲁 棒 性 ， 同 时 
丙 斯 电 启 。 001 08360 07592 09428 09917 均衡 不 可 见 性 ， 可 以 有 效 提高 知识 产权 的 保护 和 认证 。 
椒盐 噪声 0.01 0.8560 0.8953 0.9286 1.0000 
高 斯 白 噪声 20 0.8428 0.9012 0.9429 0.9944 参考 文献 : 
旋转 玫 击 网 人 [1] 胡 裕 峰 . 图 像 变换 数字 水 印 技术 研究 [D]. 杭州 : 浙江 大 学 , 2014. 
人 [2] 杨 杭 馈 ,玉玲 . 非 执 样 Contourlet 变换 的 低频 数字 水 印 [J] 计算 机 工 
图 像 变 暗 (0，0.92) 0.8371 0.8326 0.9328 0.9861 程 与 应 用 2011, 47 (16): 201-203. 
缩放 攻击 1/2 0.9407 0.8872 0.9000 0.9848 
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